Spin-fluctuations in the quarter-filled Hubbard ring : significances to LiV 2 04 
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Using the quantum Monte Carlo method, we investigate the spin dynamics of itinerant electrons 
in the one-dimensional Hubbard system. Based on the model calculation, we have studied the 
spin-fluctuations in the quarter-filled metallic Hubbard ring, which is aimed at the vanadium ring 
or chain defined along corner- sharing tetrahedra of LiV2 04, and found the dramatic changes of 
magnetic responses and spin-fluctuation characteristics with the temperature. Such results can 
explain the central findings in the recent neutron scattering experiment for LiV2 04. 



L1V2O4 is a transition metal oxide with a cubic spinel 
(or pyrochlore) structure showing many essential features 
of the heavy- fermion system like Ce compounds 0] . Its 
specific heat coefficient is the largest one observed among 
other 3d metallic systems, 7 ~ 420mJ/molK 2 . It has 
become an important issue to clarify the physical ori- 
gin of a high density of low-energy fermionic excitations 
without localized /-levels. There's also a great inter- 
est in the unusual magnetic properties of LiV2C>4 due 
to the itinerant frustrated nature. Another metallic sys- 
tem, Y(Sc)Mn 2 |^], has also been focused on as a frus- 
trated spin liquid belonging to a similar class to L1V2O4. 
Y(Sc)Mn2 exhibits many similarities; it is the geometri- 
cally frustrated magnet with no long range order, nearly 
antiferromagnetic itinerant system, and most interest- 
ingly the heavy fermion system. 

Three classes of theoretical mechanisms are most fre- 
quently referred to understand the heavy-fermion prop- 
erties in L1V2O4. One is the spin- fluctuations in the 
three-dimensional frustrated lattice as in the study of 
Y(Sc)Mn2H- Due to the magnetic frustrations, the spin 
cannot order down to low temperatures, resulting in the 
large fermionic entropy. The other is the well-known 
Kondo effect. From the band structure calculations^, |), 
it is shown that 3d t^g bands of V are crossing the Fermi 
level and, by the trigonal crystal field, split to a bit 
narrow-band half-filled A\ g singlet and a bit wide-band 
quarter-filled E g doublet. These results lead to the map- 
ping of the electronic structure into the Kondo lattice 
modelQ. The third candidate is the mechanism based 
on the one-dimensional electronic structure, where it is 
expected that the correlation effect is much enhanced, 
giving the large specific- heat coefficient. Fulde et aZ.@ 
have suggested that the large 7 coefficient results from 
excitations of Heisenberg spin 1/2 chains and rings, which 
are by the direct consequence of the frustration of corner- 
sharing tetrahedra of the vanadium lattice. Their idea 
on the formation of spin chain or ring in LiV204 dates 
back to the study on Yb 4 As3||7|], where, due to a charge 
ordering of Yb ions, the electronic structure could be in- 
terpreted as well-decoupled (at least magnetically) one- 
dimensional chains. 

In the last decades, the spin-fluctuation has been 
found to play fundamental roles in many of strongly 
correlated electron systems, especially in the high Tc 
superconductors M. The combined system of the strong 



electron correlation and the spin itineracy, where there is 
no magnetic long range order, leads to intriguing spin- 
fluctuations. Recent two inelastic neutron scattering 
experiments^, |l^| have delivered seminal informations 
on the spin dynamics in L1V2O4; (i) they have reported 
the dramatic crossover from a ferromagnetic (FM) to an 
antiferromagnetic (AFM) spin-fluctuation with the tem- 
perature T lowered. Especially, Lee et aZ.jHj have ex- 
plicitly pointed out the AFM spin-fluctuation would be 
centered around Q c — 0.64A -1 = 0.59ir/a (a is the V-V 
distance), (ii) they have found the residual relaxation rate 
for T — > and its monotonous increase at Q c with raising 
temperature, but the increasing behavior was reported 
differently from each other, i.e. Krimmel et a/.j|] have 
reported the square-root temperature behavior, whereas 
Lee et q/.|lO|] the linear behavior, and (hi) Krimmel et al. 
have also provided the momentum-transfer-dependence 
of the relaxation rate to elucidate the change of spin- 
fluctuation characteristics. 

In this paper, we bring focus on the spin dynamics of 
L1V2O4 grounded on the one-dimensional mechanism. It 
is actually clear that the Heisenberg spin 1/2 chain by 
itself cannot explain the observed experiments because 
the low energy excitation of the system should be well- 
dispersive gapless magnon. Instead, the quarter-filled 
Hubbard ring is taken as the starting point, which gives 
the itineracy to Fulde's spins, toward an understanding 
of characteristic spin- fluctuations observed in L1V2O4. 
Fujimoto|]Tl| has studied the network of quarter-filled 
Hubbard chains accounting for the hybridization between 
chains along the similar line to the Fulde's spin chain or 
ring. In the study, he has obtained quite a comparable 
size of 7 to the experimental value and introduced an- 
other energy scale T* giving the dimensional crossover; 
below T*, three-dimensional Fermi-liquid state with the 
heavy mass is realized, but above T* , one-dimensional 
characters dominate. He has also pointed out that Urano 
et al. fl2]| 's transport data can be consistent with this pic- 
ture and a characteristic low temperature scale (~ 20 
K) observed would correspond to T* . In the present in- 
vestigation, therefore, we do not consider the very low 
temperature region of T -C T* (li| . 

In an actual situation of L1V2O4, the one-dimensional 
ring (chain) is constructed on the corner-sharing tetrahe- 
dra of the vanadium network, which has been visualized 
in Refs.|6|, O. The Hubbard model is defined on the 
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one-dimensional lattice; 



i S=±l a i 

+ U ^n^riii, (1) 



where c\ a and Ci a are the creation and annihilation opera- 
tors for electrons with spin a at lattice i and rii a — c[ a Ci a . 
t is the hopping parameter, fj, the chemical potential, and 
U the on-site Coulomb correlation. The first and the last 
sites are connected by imposing the periodic boundary 
condition, i.e. it gives the "ring" geometry. 

Among many versions of quantum Monte Carlo 
(QMC) methods, it is the path integral theory of QMC 
that is better proper for a description of the itinerant 
electron systems [[5], The path integral QMC for 
the correlated electron system cooperates the Hubbard- 
Stratonovitch transformation and integrates out the elec- 
tron field. Most of all the QMC methods are based on 
the Trotter decomposition e -P(K+v) K ^-AtV e -ArK^L 

with (3 = 1/T = ArL0 where K should be the one- 
electron terms and V the electron-electron correlation 
term. The collective spin excitations probed by the 
inelastic neutron scattering are described in the time- 
dependent spin-spin correlation function 5(q, r) 



S(q,r) 



J_ \p iq.(R 4 -R 3 ) 



[ni T (T)-n < j.(r)][n3 T -n i J), 

which is corresponding to the thermodynamic two- 
particle Green's function in the imaginary time. Through 
the analytic continuation from S^q, r) under a condition 
5(q, u>) > 0[[l8|, we obtain the experimentally observable 
spectral function 5(q, u>) satisfying 
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For a numerical simulation in the study, the 24-sitc 
Hubbard ring with U/t = 4 is considered. We follow the 
basic approach for the grand canonical ensemble. For 
the quarter-filled occupation, the ensembles such that it 
can give ^/NJ2i( n tf + n u) = 0-5 should be sampled by 
adjusting the value of /i. The Trotter decomposition is 
done such that dr = 0.1 (only for a case of 1/T = 1.5, 
dr = 0.05). We have taken the averages of the dynamical 
correlation functions over 10 4 updates of all the Hubbard- 
Stratonovitch bosons on the lattice. Further, to keep the 
numerical stability, we have used the matrix factorization 
technique Jl5[] . All the energy quantities are measured in 
a unit of t and all the momentum quantities in Tt/a. 

Figure |l|(a)-(e) show the dynamical spin-spin correla- 
tion function in the quasi-elastic mode of the quarter- 
filled Hubbard ring at various temperatures. As the tem- 
perature decreases, at first, the AFM short-range order, 
induced by the AFM spin-fluctuation, with a character- 
istic wave vector q = 1/2 becomes appreciable. On the 
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FIG. 1: (a)-(e) Dynamical spin-spin correlation function in 
the quarter-filled Hubbard ring as a function of q in the quasi- 
elastic mode (lo = 0) at several temperatures, i.e. S(q,u = 0). 
(f) Ratio of S(q = 1/2, u) = 0) and S(q = Q,u> = 0). 



other hand, in the high temperature region, its character 
is found rather FM, more obvious in Fig.^. This finding 
is very consistent with the dramatic crossover of the char- 
acter of spin-fluctuation with the temperature reported 
by recent inelastic neutron scattering experiments ||, [l0|| . 
Further, the characteristic wave vector q = 1/2 is also 
moderately consistent with the staggered wave vector 
Q c = 0.59 found in the experiment p0| . It can be un- 
derstood as follows; the one-dimensional spin-spin corre- 
lation function in the itinerant limit gives rise to a log- 
arithmic singularity at 2k p (&f is the Fermi wave vec- 
tor; kp = 1/4 in the present case) [^9) , from which in 
the correlated limit incomplete magnetic moments (i.e. 
AFM spin-fluctuation) could evolve at every other site 
and lead to the AFM short-range order. In the low tem- 
perature region, S(q, 0) is rapidly suppressed for high 
momentum transfers (q > 1/2). In Fig.|l|(f), the ratio 
of quasi-elastic peaks of the magnetic scattering at q = 
and 5 = 1/2 shows roughly the change of characteristic 
spin- fluctuations governing the system around 1/T ~ 4 
or T ~ 0.25. Here we need estimate the size of t. Fulde 
et al. H have estimated the exchange integral of J in the 
Heisenberg spin 1/2 chain on the frustrated lattice as 
about 3 meVQ. Noting J ~ 0(t 2 /U), we then have ap- 
proximately t ~ 10 meV. Now the crossover temperature 
is T ~ 0.25 - 30 K, close to Krimmel et aVs 40 K§. 

To scrutinize the changes of magnetic responses more, 
we provide, in Fig.|2|, the static susceptibility Xq at <1 = 
and q = 1/2 with T. The static susceptibility x q can 
be directly evaluated from 5(q, r) obtained by the QMC 
calculation 



X q = / rfr5(q,r). 
Jo 



(4) 
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FIG. 2: Crossover from FM to AFM spin-fluctuation with T. 
The static spin susceptibility at q — and q = 1/2; the solid 
line for q = 1/2 is the AFM Curie- Weiss susceptibility and 
the dotted line for q — is a guide for the eye. 



It is shown in the figure that \q 101 <Z = 1/2 nicely fol- 
lows the AFM Curie- Weiss susceptibility, oc 1/(T + 9). 
Such Curie- Weiss behavior has been already observed 
in the experiment ]lC| , where 9 for the best fitting was 
estimated as 7.5 K. Our calculation gives the similar 
value, 9 ~ 0.11 ~ 13.2 K. Let us note it means that 
the quarter-filled one-dimensional Hubbard model could 
allow a formation of magnetic moments like Kondo lat- 
tice model in the high temperature (2Cj, where the local- 
ized level manifests the Curie- Weiss susceptibility. Un- 
stable magnetic moments produce the AFM short-range 
order centered around w/(2a) (i.e. q = 1/2), where 2a 
is the distance between neighboring moments. It is con- 
sistent with the temperature behavior of the magnetic 
relaxation rate T q at q = 1/2 in Fig.||. Fig.^ also shows 
the crossover of magnetic characters from FM to AFM 
around T ~ 0.25 - 30 K, consistent with Figjljf). 

Another important quantity is the magnetic relaxation 
rate T q usually defined by the simple ansatz for the dy- 



namic susceptibility 



wr q xq 
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In the study, T q is evaluated by taking u> — > in Eq. (^|) , 

<?(q,0) = ^, 
1 q 

where S(q, 0) and Xq are already given in Figs.[l] and|2|, re- 
spectively. But it should be noted that, because the unit 
of S(q,0) is arbitrary pl|, T q would be obtained only up 
to a constant. That is, we note the true relaxation rate 
should be nT q (rj is a nonzero constant). The results of 
T q are provided in Fig||. The upper panel of Fig|j| shows 
that the spin magnetic relaxation rate r g at q — 1/2 in- 
creases linearly with temperature for a rather wide tem- 
perature range to ~ 0.7 ~ 80 K. The linear increasing 
behavior for such a wide T range (to ~ 80 K) has been 
ascertained in the experiment by Lee et al. (!□] , where its 
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FIG. 3: Magnetic relaxation rate F q . Upper panel: the tem- 
perature dependence of r q at q — 1/2. Lower panel: the 
q-dependence of T q at high and low temperatures. 



increasing rate is found 0.46. Linear T behavior is actu- 
ally unusual in the /-electron heavy fermion system. It 
is however a bit common in frustrated metal oxides, nor- 
mally related with unstable local moments. In the Kondo 
system, one most usually has r g ~ T° + bT 1 / 2 , where 
T q ~ Tk (Kondo temperature) Q . Earlier, through 
mapping into the Kondo model, Anisimov et al.^ have 
estimated Tk ~ 550 K for the single-site case, but argued 
that another characteristic energy scale T co ^(^ 25 — 40 K, 
comparable to T°) would replace Tk in the dense Kondo 
lattice. The increasing rate of the present result shown 
in Fig.|| is estimated as 0.32 and at T = is ob- 
tained about 0.1 meV. Comparing the increasing rates, 
we find an unknown constant 77 should be about 1.44 and 
the residual relaxation rate T° be 0.14 meV. This value 
is smaller than experimental findings, i.e. Krimmel et 
al.^ have reported 0.5 meV and Lee et q/.p0| 1.4 meV. 
The lower panel of Fig.|| gives the q dependences of T q , 
whose dependences also agree with an observation of the 
experiment^]. For 1/T = 1.5 (high temperature), T q 
shows a linear q dependence, which is actually expected 
in the spin-fluctuation theories of weak FM metals [p2[. 
Therefore, the behavior of linear q dependence is con- 
sistent with our argument that the system should be a 
metal with weak FM spin-fluctuations at high tempera- 
tures (> 0.25), being associated with Figs.|l| and 0. On 
the other hand, at the low temperature (1/T = 10) 
is almost constant for small q's (q < 1/2), but rapidly 
increases for high q's (q > 1/2). The constant T q with q 
at low T (not anticipated by the simple Fermi-liquid the- 
ory) was found in the experiment, but a rapid increase for 
high g's was not, which instead may be attributed to an- 
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other subtle feature of one-dimensional Hubbard model. 
Recently, it has been found that the nonlinear coupling 
between spin and charge in the one-dimensional Hubbard 
model would lead to coupled collective excitations other 
than spin-fluctuations (or magnons) in S(q, u>) especially 
for high g's[|l^|. Those may serve as additional decaying 
channels for spin fluctuations. It is noted that such en- 
hanced T q is directly connected with the diminution of 
S(q, 0) at low T for q > 1/2 in Fig|[ 

In summary, we have discussed the recent inelastic 
neutron scattering experiments for L1V2O4 based on the 
QMC study of spin-fluctuations in the quarter-filled Hub- 
bard ring. In the study, neutron scattering cross sections, 
static spin susceptibilities, and spin relaxation rates have 
been evaluated with temperatures and momentum trans- 
fers. They are found quite consistent with the experi- 
ment qualitatively, or semi-quantitatively. Particularly, 
it is appealing that the AFM short-range correlation de- 
velops due to a formation of unstable magnetic moments 
as T decreases in the quarter-filled Hubbard ring, which 
explains the AFM spin-fluctuation around Q c = 0.59 ob- 
served in L1V2O4. Our finding that a single quarter- 
filled Hubbard ring can explain the neutron scattering 
experiment could be consistent with a case of decou- 
pled chains in Yb4As3 unless we think of the very low 
temperature regime. However, it is a difference that the 
one-dimensional electronic structure is expected from the 
geometrical frustration in LiV204jn}. We would like to 
remark that at least a few experimental findings cannot 



be explained by the Kondo lattice model, but by the one- 
dimensional quarter-filled metallic model; (i) features of 
nearly-FM metal at high temperatures, (ii) the linear T 
dependence of T q for a wide T range, and (iii) the con- 
stant T g with q at low temperatures (actually decreas- 
ing behaviors in CeCue(2^]). Further, the recent nuclear 
magnetic resonance (NMR) study for L1V2O4 under high 
pressure has reported an opposite behavior of T\ (spin- 
lattice relaxation time) to that of Ce compounds ]24|. It 
is also worth stressing that magnetic responses (such as 
T dependence of x) °f L1V2O4 differ qualitatively from 
Y(Sc)Mn2^]. Therefore, the present results can be one 
of evidences along with other studies || Q that LIV2O4 
comprises one-dimensional chains or rings and behaves 
like the one-dimensional system. Finally, it should be 
noted that the conclusion casts another important prob- 
lem to us. It is well known that, in one-dimensional 
metallic system, some extreme realization of correlation 
effects like the spin-charge separation occurs, which has 
been actually observed in SrCuC>2 by the photoemission 
spectroscopy Thus it must be fascinating to search 
for the spin-charge separation in LiV204 by the photoe- 
mission spectroscopy. 
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